<!-- saved from url=(0014)about:internet -->
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8" />
<title>RTW Report - ode78_.c</title>
<script language="JavaScript" type="text/javascript" src="rtwreport_utils.js"></script>
</head>
<body bgcolor="#eeeeee" text="#1122aa" onload="try {if (top) { if (top.rtwFileOnLoad) top.rtwFileOnLoad(window.document); else local_onload();}} catch(err) {};">
<p>
<table border="0" cellspacing="0" cellpadding="6" width="100%%" height="100%%"><tr><td width="100%%" valign="top" bgcolor="#ffffff">
<h4>File: <a href="../ode78_.c" target="rtwreport_document_frame" id="linkToText_plain">ode78_.c</a></h4>
<pre id="RTWcode">
<a name="1"><font color="#888888 "><i>    1   </i></font></a><font color="#117755"><i>/*</i></font>
<a name="2"><font color="#888888 "><i>    2   </i></font></a><font color="#117755"><I> * ode78_.c</i></font>
<a name="3"><font color="#888888 "><i>    3   </i></font></a><font color="#117755"><I> *</i></font>
<a name="4"><font color="#888888 "><i>    4   </i></font></a><font color="#117755"><I> * Code generation for function 'ode78_'</i></font>
<a name="5"><font color="#888888 "><i>    5   </i></font></a><font color="#117755"><I> *</i></font>
<a name="6"><font color="#888888 "><i>    6   </i></font></a><font color="#117755"><I> * C source code generated on: Mon Jun 22 23:07:02 2015</i></font>
<a name="7"><font color="#888888 "><i>    7   </i></font></a><font color="#117755"><I> *</i></font>
<a name="8"><font color="#888888 "><i>    8   </i></font></a><font color="#117755"><I> */</i></font>
<a name="9"><font color="#888888 "><i>    9   </i></font></a>
<a name="10"><font color="#888888 "><i>   10   </i></font></a><font color="#117755"><i>/* Include files */</i></font>
<a name="11"><font color="#888888 "><i>   11   </i></font></a><font color="#992211">#</font><font color="#992211"><b>include</b></font> <font color="#992211">&quot;rt_nonfinite.h&quot;</font>
<a name="12"><font color="#888888 "><i>   12   </i></font></a><font color="#992211">#</font><font color="#992211"><b>include</b></font> <font color="#992211">&quot;ode78_.h&quot;</font>
<a name="13"><font color="#888888 "><i>   13   </i></font></a><font color="#992211">#</font><font color="#992211"><b>include</b></font> <font color="#992211">&quot;fy_.h&quot;</font>
<a name="14"><font color="#888888 "><i>   14   </i></font></a><font color="#992211">#</font><font color="#992211"><b>include</b></font> <font color="#992211">&quot;ode78__data.h&quot;</font>
<a name="15"><font color="#888888 "><i>   15   </i></font></a>
<a name="16"><font color="#888888 "><i>   16   </i></font></a><font color="#117755"><i>/* Variable Definitions */</i></font>
<a name="17"><font color="#888888 "><i>   17   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_emlrtRSI">emlrtRSI</a> = <b>{</b> 87, <font color="#1122ff">&quot;ode78_&quot;</font>,
<a name="18"><font color="#888888 "><i>   18   </i></font></a>  <font color="#1122ff">&quot;C:/\xe4\xba\x91\xe7\x9b\x98\xe5\x90\x8c\xe6\xad\xa5\xe6\x96\x87\xe4\xbb\xb6\xe5\xa4?/GTOC8/Lowthrust_mee_v0/ode78_.m&quot;</font>
<a name="19"><font color="#888888 "><i>   19   </i></font></a><b>}</b>;
<a name="20"><font color="#888888 "><i>   20   </i></font></a>
<a name="21"><font color="#888888 "><i>   21   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_b_emlrtRSI">b_emlrtRSI</a> = <b>{</b> 89, <font color="#1122ff">&quot;ode78_&quot;</font>,
<a name="22"><font color="#888888 "><i>   22   </i></font></a>  <font color="#1122ff">&quot;C:/\xe4\xba\x91\xe7\x9b\x98\xe5\x90\x8c\xe6\xad\xa5\xe6\x96\x87\xe4\xbb\xb6\xe5\xa4?/GTOC8/Lowthrust_mee_v0/ode78_.m&quot;</font>
<a name="23"><font color="#888888 "><i>   23   </i></font></a><b>}</b>;
<a name="24"><font color="#888888 "><i>   24   </i></font></a>
<a name="25"><font color="#888888 "><i>   25   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_c_emlrtRSI">c_emlrtRSI</a> = <b>{</b> 129, <font color="#1122ff">&quot;ode78_&quot;</font>,
<a name="26"><font color="#888888 "><i>   26   </i></font></a>  <font color="#1122ff">&quot;C:/\xe4\xba\x91\xe7\x9b\x98\xe5\x90\x8c\xe6\xad\xa5\xe6\x96\x87\xe4\xbb\xb6\xe5\xa4?/GTOC8/Lowthrust_mee_v0/ode78_.m&quot;</font>
<a name="27"><font color="#888888 "><i>   27   </i></font></a><b>}</b>;
<a name="28"><font color="#888888 "><i>   28   </i></font></a>
<a name="29"><font color="#888888 "><i>   29   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_i_emlrtRSI">i_emlrtRSI</a> = <b>{</b> 37, <font color="#1122ff">&quot;mpower&quot;</font>,
<a name="30"><font color="#888888 "><i>   30   </i></font></a>  <font color="#1122ff">&quot;C:/Program Files/MATLAB/R2013a/toolbox/eml/lib/matlab/ops/mpower.m&quot;</font> <b>}</b>;
<a name="31"><font color="#888888 "><i>   31   </i></font></a>
<a name="32"><font color="#888888 "><i>   32   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_j_emlrtRSI">j_emlrtRSI</a> = <b>{</b> 42, <font color="#1122ff">&quot;power&quot;</font>,
<a name="33"><font color="#888888 "><i>   33   </i></font></a>  <font color="#1122ff">&quot;C:/Program Files/MATLAB/R2013a/toolbox/eml/lib/matlab/ops/power.m&quot;</font> <b>}</b>;
<a name="34"><font color="#888888 "><i>   34   </i></font></a>
<a name="35"><font color="#888888 "><i>   35   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_k_emlrtRSI">k_emlrtRSI</a> = <b>{</b> 56, <font color="#1122ff">&quot;power&quot;</font>,
<a name="36"><font color="#888888 "><i>   36   </i></font></a>  <font color="#1122ff">&quot;C:/Program Files/MATLAB/R2013a/toolbox/eml/lib/matlab/ops/power.m&quot;</font> <b>}</b>;
<a name="37"><font color="#888888 "><i>   37   </i></font></a>
<a name="38"><font color="#888888 "><i>   38   </i></font></a><font color="#117755"><i>/* Function Declarations */</i></font>
<a name="39"><font color="#888888 "><i>   39   </i></font></a><font color="#112266"><b>static</b></font> <font color="#112266"><b>void</b></font> <a href="#fcn_b_eml_error" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_b_eml_error');" target="rtwreport_document_frame"><font color="#1122aa">b_eml_error</font></a>(<font color="#112266"><b>void</b></font>);
<a name="40"><font color="#888888 "><i>   40   </i></font></a>
<a name="41"><font color="#888888 "><i>   41   </i></font></a><font color="#117755"><i>/* Function Definitions */</i></font>
<a name="42"><font color="#888888 "><i>   42   </i></font></a><font color="#112266"><b>static</b></font> <font color="#112266"><b>void</b></font> <a name="fcn_b_eml_error">b_eml_error</a>(<font color="#112266"><b>void</b></font>)
<a name="43"><font color="#888888 "><i>   43   </i></font></a><b>{</b>
<a name="44"><font color="#888888 "><i>   44   </i></font></a>  emlrtPushRtStackR2012b(&amp;<a href="ode78__data_c.html#var_h_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_h_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">h_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="45"><font color="#888888 "><i>   45   </i></font></a>  emlrtErrorWithMessageIdR2012b(<a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>, &amp;<a href="ode78__data_c.html#var_emlrtRTEI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRTEI');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRTEI</font></a>,
<a name="46"><font color="#888888 "><i>   46   </i></font></a>    <font color="#1122ff">&quot;Coder:toolbox:power_domainError&quot;</font>, 0);
<a name="47"><font color="#888888 "><i>   47   </i></font></a>  emlrtPopRtStackR2012b(&amp;<a href="ode78__data_c.html#var_h_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_h_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">h_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="48"><font color="#888888 "><i>   48   </i></font></a><b>}</b>
<a name="49"><font color="#888888 "><i>   49   </i></font></a>
<a name="50"><font color="#888888 "><i>   50   </i></font></a><font color="#112266"><b>void</b></font> <a name="fcn_ode78_">ode78_</a>(<font color="#112266"><b>const</b></font> real_T tspan[2], <font color="#112266"><b>const</b></font> real_T x0[14], <font color="#112266"><b>const</b></font> <a href="ode78__types_h.html#type_struct_T" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'type_struct_T');" target="rtwreport_document_frame"><font color="#1122aa">struct_T</font></a> *auxdata,
<a name="51"><font color="#888888 "><i>   51   </i></font></a>            real_T *tout, real_T xout[14], real_T *ErrorFlag)
<a name="52"><font color="#888888 "><i>   52   </i></font></a><b>{</b>
<a name="53"><font color="#888888 "><i>   53   </i></font></a>  real_T hmax;
<a name="54"><font color="#888888 "><i>   54   </i></font></a>  real_T h;
<a name="55"><font color="#888888 "><i>   55   </i></font></a>  real_T t;
<a name="56"><font color="#888888 "><i>   56   </i></font></a>  int32_T fk;
<a name="57"><font color="#888888 "><i>   57   </i></font></a>  real_T x[14];
<a name="58"><font color="#888888 "><i>   58   </i></font></a>  real_T f[182];
<a name="59"><font color="#888888 "><i>   59   </i></font></a>  int32_T i;
<a name="60"><font color="#888888 "><i>   60   </i></font></a>  int32_T i0;
<a name="61"><font color="#888888 "><i>   61   </i></font></a>  boolean_T exitg1;
<a name="62"><font color="#888888 "><i>   62   </i></font></a>  real_T b[13];
<a name="63"><font color="#888888 "><i>   63   </i></font></a>  <font color="#112266"><b>static</b></font> <font color="#112266"><b>const</b></font> real_T a_i_j[156] = <b>{</b> 0.055555555555555552, 0.0, 0.0, 0.0, 0.0,
<a name="64"><font color="#888888 "><i>   64   </i></font></a>    0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.020833333333333332, 0.0625, 0.0,
<a name="65"><font color="#888888 "><i>   65   </i></font></a>    0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03125, 0.0, 0.09375, 0.0,
<a name="66"><font color="#888888 "><i>   66   </i></font></a>    0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.3125, 0.0, -1.171875,
<a name="67"><font color="#888888 "><i>   67   </i></font></a>    1.171875, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0375, 0.0, 0.0,
<a name="68"><font color="#888888 "><i>   68   </i></font></a>    0.1875, 0.15, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.047910137111111112,
<a name="69"><font color="#888888 "><i>   69   </i></font></a>    0.0, 0.0, 0.11224871277777777, -0.025505673777777779, 0.012846823888888888,
<a name="70"><font color="#888888 "><i>   70   </i></font></a>    0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.016917989787292281, 0.0, 0.0,
<a name="71"><font color="#888888 "><i>   71   </i></font></a>    0.3878482784860432, 0.035977369851500331, 0.19697021421566607,
<a name="72"><font color="#888888 "><i>   72   </i></font></a>    -0.17271385234050185, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0690957533591923, 0.0,
<a name="73"><font color="#888888 "><i>   73   </i></font></a>    0.0, -0.63424797672885413, -0.16119757522460407, 0.13865030945882525,
<a name="74"><font color="#888888 "><i>   74   </i></font></a>    0.94092861403575623, 0.21163632648194397, 0.0, 0.0, 0.0, 0.0, 0.0,
<a name="75"><font color="#888888 "><i>   75   </i></font></a>    0.18355699683904539, 0.0, 0.0, -2.4687680843155926, -0.29128688781630047,
<a name="76"><font color="#888888 "><i>   76   </i></font></a>    -0.026473020233117376, 2.8478387641928005, 0.28138733146984979,
<a name="77"><font color="#888888 "><i>   77   </i></font></a>    0.12374489986331466, 0.0, 0.0, 0.0, 0.0, -1.2154248173958881, 0.0, 0.0,
<a name="78"><font color="#888888 "><i>   78   </i></font></a>    16.672608665945774, 0.915741828416818, -6.0566058043574706,
<a name="79"><font color="#888888 "><i>   79   </i></font></a>    -16.00357359415618, 14.849303086297663, -13.371575735289849,
<a name="80"><font color="#888888 "><i>   80   </i></font></a>    5.134182648179638, 0.0, 0.0, 0.0, 0.25886091643826425, 0.0, 0.0,
<a name="81"><font color="#888888 "><i>   81   </i></font></a>    -4.7744857854892047, -0.43509301377703252, -3.0494833320722416,
<a name="82"><font color="#888888 "><i>   82   </i></font></a>    5.5779200399360995, 6.15583158986104, -5.0621045867369387, 2.193926173180679,
<a name="83"><font color="#888888 "><i>   83   </i></font></a>    0.13462799865933495, 0.0, 0.0, 0.82242759962650747, 0.0, 0.0,
<a name="84"><font color="#888888 "><i>   84   </i></font></a>    -11.658673257277664, -0.75762211669093615, 0.71397358815958156,
<a name="85"><font color="#888888 "><i>   85   </i></font></a>    12.075774986890057, -2.1276591139204029, 1.9901662070489554,
<a name="86"><font color="#888888 "><i>   86   </i></font></a>    -0.23428647154404028, 0.17589857770794226, 0.0, 0.0 <b>}</b>;
<a name="87"><font color="#888888 "><i>   87   </i></font></a>
<a name="88"><font color="#888888 "><i>   88   </i></font></a>  real_T b_x[14];
<a name="89"><font color="#888888 "><i>   89   </i></font></a>  real_T scale;
<a name="90"><font color="#888888 "><i>   90   </i></font></a>  int32_T i1;
<a name="91"><font color="#888888 "><i>   91   </i></font></a>  real_T sol2[14];
<a name="92"><font color="#888888 "><i>   92   </i></font></a>  <font color="#112266"><b>static</b></font> <font color="#112266"><b>const</b></font> real_T b_b[13] = <b>{</b> 0.041747491141530244, 0.0, 0.0, 0.0, 0.0,
<a name="93"><font color="#888888 "><i>   93   </i></font></a>    -0.055452328611239311, 0.23931280720118009, 0.703510669403443,
<a name="94"><font color="#888888 "><i>   94   </i></font></a>    -0.75975961381446089, 0.6605630309222863, 0.15818748251012332,
<a name="95"><font color="#888888 "><i>   95   </i></font></a>    -0.23810953875286281, 0.25 <b>}</b>;
<a name="96"><font color="#888888 "><i>   96   </i></font></a>
<a name="97"><font color="#888888 "><i>   97   </i></font></a>  <font color="#112266"><b>static</b></font> <font color="#112266"><b>const</b></font> real_T c_b[13] = <b>{</b> 0.0295532136763535, 0.0, 0.0, 0.0, 0.0,
<a name="98"><font color="#888888 "><i>   98   </i></font></a>    -0.828606276487797, 0.31124090005111832, 2.4673451905998869,
<a name="99"><font color="#888888 "><i>   99   </i></font></a>    -2.5469416518419088, 1.4435485836767752, 0.079415595881127288,
<a name="100"><font color="#888888 "><i>  100   </i></font></a>    0.044444444444444446, 0.0 <b>}</b>;
<a name="101"><font color="#888888 "><i>  101   </i></font></a>
<a name="102"><font color="#888888 "><i>  102   </i></font></a>  real_T error_1;
<a name="103"><font color="#888888 "><i>  103   </i></font></a>  real_T absx;
<a name="104"><font color="#888888 "><i>  104   </i></font></a>  real_T b_t;
<a name="105"><font color="#888888 "><i>  105   </i></font></a>  boolean_T exitg2;
<a name="106"><font color="#888888 "><i>  106   </i></font></a>
<a name="107"><font color="#888888 "><i>  107   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------ */</i></font>
<a name="108"><font color="#888888 "><i>  108   </i></font></a>  <font color="#117755"><i>/*  rk78 integration */</i></font>
<a name="109"><font color="#888888 "><i>  109   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="110"><font color="#888888 "><i>  110   </i></font></a>  <font color="#117755"><i>/*  [tout, xout, ErrorFlag] = ode78_chen_(tspan , x0 , auxdata) */</i></font>
<a name="111"><font color="#888888 "><i>  111   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="112"><font color="#888888 "><i>  112   </i></font></a>  <font color="#117755"><i>/*  Input argumuents: */</i></font>
<a name="113"><font color="#888888 "><i>  113   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------------- */</i></font>
<a name="114"><font color="#888888 "><i>  114   </i></font></a>  <font color="#117755"><i>/*  tspan            [1x2 double]           time array                                  [ TU ]         */</i></font>
<a name="115"><font color="#888888 "><i>  115   </i></font></a>  <font color="#117755"><i>/*  x0                 [14x1 double]         state and costate variables     [  -  ] */</i></font>
<a name="116"><font color="#888888 "><i>  116   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="117"><font color="#888888 "><i>  117   </i></font></a>  <font color="#117755"><i>/*  Output argumuents: */</i></font>
<a name="118"><font color="#888888 "><i>  118   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------------- */</i></font>
<a name="119"><font color="#888888 "><i>  119   </i></font></a>  <font color="#117755"><i>/*  tout              [1x1 double]            position vector          [  -  ] */</i></font>
<a name="120"><font color="#888888 "><i>  120   </i></font></a>  <font color="#117755"><i>/*  xout             [14x1 double]          velocity vector           [  -  ] */</i></font>
<a name="121"><font color="#888888 "><i>  121   </i></font></a>  <font color="#117755"><i>/*  ErrorFlag   [1x1 double]             error flag                   [  -  ] */</i></font>
<a name="122"><font color="#888888 "><i>  122   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="123"><font color="#888888 "><i>  123   </i></font></a>  <font color="#117755"><i>/*  External functions called: */</i></font>
<a name="124"><font color="#888888 "><i>  124   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------------- */</i></font>
<a name="125"><font color="#888888 "><i>  125   </i></font></a>  <font color="#117755"><i>/*  fy_.m */</i></font>
<a name="126"><font color="#888888 "><i>  126   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="127"><font color="#888888 "><i>  127   </i></font></a>  <font color="#117755"><i>/*  Copyright (C) 2015/06/22 by Chen Zhang,  */</i></font>
<a name="128"><font color="#888888 "><i>  128   </i></font></a>  <font color="#117755"><i>/*  BeiHang University, chenzhang.buaa@gmail.com */</i></font>
<a name="129"><font color="#888888 "><i>  129   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------ */</i></font>
<a name="130"><font color="#888888 "><i>  130   </i></font></a>  <font color="#117755"><i>/*  The coefficients of rk78 */</i></font>
<a name="131"><font color="#888888 "><i>  131   </i></font></a>  <font color="#117755"><i>/*  power for step control */</i></font>
<a name="132"><font color="#888888 "><i>  132   </i></font></a>  <font color="#117755"><i>/*  Maximum step size */</i></font>
<a name="133"><font color="#888888 "><i>  133   </i></font></a>  hmax = (tspan[1] - tspan[0]) / 2.5;
<a name="134"><font color="#888888 "><i>  134   </i></font></a>
<a name="135"><font color="#888888 "><i>  135   </i></font></a>  <font color="#117755"><i>/*  Initial step size */</i></font>
<a name="136"><font color="#888888 "><i>  136   </i></font></a>  h = hmax;
<a name="137"><font color="#888888 "><i>  137   </i></font></a>  <font color="#112266"><b>if</b></font> (hmax &gt; 0.5) <b>{</b>
<a name="138"><font color="#888888 "><i>  138   </i></font></a>    h = 0.5;
<a name="139"><font color="#888888 "><i>  139   </i></font></a>  <b>}</b>
<a name="140"><font color="#888888 "><i>  140   </i></font></a>
<a name="141"><font color="#888888 "><i>  141   </i></font></a>  <font color="#112266"><b>if</b></font> (h &gt; hmax) <b>{</b>
<a name="142"><font color="#888888 "><i>  142   </i></font></a>    h = hmax;
<a name="143"><font color="#888888 "><i>  143   </i></font></a>  <b>}</b>
<a name="144"><font color="#888888 "><i>  144   </i></font></a>
<a name="145"><font color="#888888 "><i>  145   </i></font></a>  <font color="#117755"><i>/*  Relative error tolerance (applies to all components of the solution vector) */</i></font>
<a name="146"><font color="#888888 "><i>  146   </i></font></a>  t = tspan[0];
<a name="147"><font color="#888888 "><i>  147   </i></font></a>
<a name="148"><font color="#888888 "><i>  148   </i></font></a>  <font color="#117755"><i>/*  Minimum step size */</i></font>
<a name="149"><font color="#888888 "><i>  149   </i></font></a>  <font color="#117755"><i>/*  Initialization */</i></font>
<a name="150"><font color="#888888 "><i>  150   </i></font></a>  <font color="#117755"><i>/*  constant for step rejection */</i></font>
<a name="151"><font color="#888888 "><i>  151   </i></font></a>  <font color="#117755"><i>/*  counter for function evaluations */</i></font>
<a name="152"><font color="#888888 "><i>  152   </i></font></a>  fk = 0;
<a name="153"><font color="#888888 "><i>  153   </i></font></a>  *ErrorFlag = 0.0;
<a name="154"><font color="#888888 "><i>  154   </i></font></a>
<a name="155"><font color="#888888 "><i>  155   </i></font></a>  <font color="#117755"><i>/*  array f for RHS calculation */</i></font>
<a name="156"><font color="#888888 "><i>  156   </i></font></a>  *tout = tspan[0];
<a name="157"><font color="#888888 "><i>  157   </i></font></a>  <font color="#112266"><b>for</b></font> (i = 0; i &lt; 14; i++) <b>{</b>
<a name="158"><font color="#888888 "><i>  158   </i></font></a>    x[i] = x0[i];
<a name="159"><font color="#888888 "><i>  159   </i></font></a>
<a name="160"><font color="#888888 "><i>  160   </i></font></a>    <font color="#117755"><i>/*  start point */</i></font>
<a name="161"><font color="#888888 "><i>  161   </i></font></a>    <font color="#112266"><b>for</b></font> (i0 = 0; i0 &lt; 13; i0++) <b>{</b>
<a name="162"><font color="#888888 "><i>  162   </i></font></a>      f[i + 14 * i0] = x0[i] * 0.0;
<a name="163"><font color="#888888 "><i>  163   </i></font></a>    <b>}</b>
<a name="164"><font color="#888888 "><i>  164   </i></font></a>
<a name="165"><font color="#888888 "><i>  165   </i></font></a>    xout[i] = x0[i];
<a name="166"><font color="#888888 "><i>  166   </i></font></a>  <b>}</b>
<a name="167"><font color="#888888 "><i>  167   </i></font></a>
<a name="168"><font color="#888888 "><i>  168   </i></font></a>  <font color="#117755"><i>/*  tau = tol * max(norm(x,'inf'), 1) ;  % accuracy */</i></font>
<a name="169"><font color="#888888 "><i>  169   </i></font></a>  <font color="#117755"><i>/*  The main loop */</i></font>
<a name="170"><font color="#888888 "><i>  170   </i></font></a>  <font color="#117755"><i>/*  count = 1; */</i></font>
<a name="171"><font color="#888888 "><i>  171   </i></font></a>  exitg1 = FALSE;
<a name="172"><font color="#888888 "><i>  172   </i></font></a>  <font color="#112266"><b>while</b></font> ((exitg1 == FALSE) &amp;&amp; ((t &lt; tspan[1]) &amp;&amp; (h &gt;= 3.5527136788005009E-15 *
<a name="173"><font color="#888888 "><i>  173   </i></font></a>           muDoubleScalarAbs(tspan[0])))) <b>{</b>
<a name="174"><font color="#888888 "><i>  174   </i></font></a>    <font color="#112266"><b>if</b></font> (t + h &gt; tspan[1]) <b>{</b>
<a name="175"><font color="#888888 "><i>  175   </i></font></a>      h = tspan[1] - t;
<a name="176"><font color="#888888 "><i>  176   </i></font></a>    <b>}</b>
<a name="177"><font color="#888888 "><i>  177   </i></font></a>
<a name="178"><font color="#888888 "><i>  178   </i></font></a>    <font color="#117755"><i>/*  Compute the RHS for step of method */</i></font>
<a name="179"><font color="#888888 "><i>  179   </i></font></a>    emlrtPushRtStackR2012b(&amp;<a href="#var_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="180"><font color="#888888 "><i>  180   </i></font></a>    <a href="fy__c.html#fcn_fy_" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_fy_');" target="rtwreport_document_frame"><font color="#1122aa">fy_</font></a>(x, auxdata-&gt;mu, auxdata-&gt;Tmax, auxdata-&gt;c, *(real_T (*)[14])&amp;f[0]);
<a name="181"><font color="#888888 "><i>  181   </i></font></a>    emlrtPopRtStackR2012b(&amp;<a href="#var_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="182"><font color="#888888 "><i>  182   </i></font></a>    <font color="#112266"><b>for</b></font> (i = 0; i &lt; 12; i++) <b>{</b>
<a name="183"><font color="#888888 "><i>  183   </i></font></a>      emlrtPushRtStackR2012b(&amp;<a href="#var_b_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_b_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">b_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="184"><font color="#888888 "><i>  184   </i></font></a>      memcpy(&amp;b[0], &amp;a_i_j[13 * i], 13U * <font color="#112266"><b>sizeof</b></font>(real_T));
<a name="185"><font color="#888888 "><i>  185   </i></font></a>      <font color="#112266"><b>for</b></font> (i0 = 0; i0 &lt; 14; i0++) <b>{</b>
<a name="186"><font color="#888888 "><i>  186   </i></font></a>        scale = 0.0;
<a name="187"><font color="#888888 "><i>  187   </i></font></a>        <font color="#112266"><b>for</b></font> (i1 = 0; i1 &lt; 13; i1++) <b>{</b>
<a name="188"><font color="#888888 "><i>  188   </i></font></a>          scale += h * f[i0 + 14 * i1] * b[i1];
<a name="189"><font color="#888888 "><i>  189   </i></font></a>        <b>}</b>
<a name="190"><font color="#888888 "><i>  190   </i></font></a>
<a name="191"><font color="#888888 "><i>  191   </i></font></a>        b_x[i0] = x[i0] + scale;
<a name="192"><font color="#888888 "><i>  192   </i></font></a>      <b>}</b>
<a name="193"><font color="#888888 "><i>  193   </i></font></a>
<a name="194"><font color="#888888 "><i>  194   </i></font></a>      <a href="fy__c.html#fcn_fy_" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_fy_');" target="rtwreport_document_frame"><font color="#1122aa">fy_</font></a>(b_x, auxdata-&gt;mu, auxdata-&gt;Tmax, auxdata-&gt;c, *(real_T (*)[14])&amp;f[14 *
<a name="195"><font color="#888888 "><i>  195   </i></font></a>          (1 + i)]);
<a name="196"><font color="#888888 "><i>  196   </i></font></a>      emlrtPopRtStackR2012b(&amp;<a href="#var_b_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_b_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">b_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="197"><font color="#888888 "><i>  197   </i></font></a>      emlrtBreakCheckFastR2012b(<a href="ode78__data_c.html#var_emlrtBreakCheckR2012bFlagVar" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtBreakCheckR2012bFlagVar');" target="rtwreport_document_frame"><font color="#1122aa">emlrtBreakCheckR2012bFlagVar</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="198"><font color="#888888 "><i>  198   </i></font></a>    <b>}</b>
<a name="199"><font color="#888888 "><i>  199   </i></font></a>
<a name="200"><font color="#888888 "><i>  200   </i></font></a>    <font color="#117755"><i>/* %%%%%%%%% Forced exit condition %%%%%%%%%% */</i></font>
<a name="201"><font color="#888888 "><i>  201   </i></font></a>    fk += 13;
<a name="202"><font color="#888888 "><i>  202   </i></font></a>    <font color="#112266"><b>if</b></font> (fk &gt; 50000) <b>{</b>
<a name="203"><font color="#888888 "><i>  203   </i></font></a>      *ErrorFlag = 1.0;
<a name="204"><font color="#888888 "><i>  204   </i></font></a>      exitg1 = TRUE;
<a name="205"><font color="#888888 "><i>  205   </i></font></a>    <b>}</b> <font color="#112266"><b>else</b></font> <b>{</b>
<a name="206"><font color="#888888 "><i>  206   </i></font></a>      <font color="#117755"><i>/* %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% */</i></font>
<a name="207"><font color="#888888 "><i>  207   </i></font></a>      <font color="#117755"><i>/*  Two solution */</i></font>
<a name="208"><font color="#888888 "><i>  208   </i></font></a>      <font color="#112266"><b>for</b></font> (i0 = 0; i0 &lt; 14; i0++) <b>{</b>
<a name="209"><font color="#888888 "><i>  209   </i></font></a>        scale = 0.0;
<a name="210"><font color="#888888 "><i>  210   </i></font></a>        <font color="#112266"><b>for</b></font> (i1 = 0; i1 &lt; 13; i1++) <b>{</b>
<a name="211"><font color="#888888 "><i>  211   </i></font></a>          scale += h * f[i0 + 14 * i1] * b_b[i1];
<a name="212"><font color="#888888 "><i>  212   </i></font></a>        <b>}</b>
<a name="213"><font color="#888888 "><i>  213   </i></font></a>
<a name="214"><font color="#888888 "><i>  214   </i></font></a>        sol2[i0] = x[i0] + scale;
<a name="215"><font color="#888888 "><i>  215   </i></font></a>      <b>}</b>
<a name="216"><font color="#888888 "><i>  216   </i></font></a>
<a name="217"><font color="#888888 "><i>  217   </i></font></a>      <font color="#117755"><i>/*  8-th order solution */</i></font>
<a name="218"><font color="#888888 "><i>  218   </i></font></a>      <font color="#117755"><i>/*  7-th order solution */</i></font>
<a name="219"><font color="#888888 "><i>  219   </i></font></a>      <font color="#117755"><i>/*  Truncation error */</i></font>
<a name="220"><font color="#888888 "><i>  220   </i></font></a>      <font color="#112266"><b>for</b></font> (i0 = 0; i0 &lt; 14; i0++) <b>{</b>
<a name="221"><font color="#888888 "><i>  221   </i></font></a>        scale = 0.0;
<a name="222"><font color="#888888 "><i>  222   </i></font></a>        <font color="#112266"><b>for</b></font> (i1 = 0; i1 &lt; 13; i1++) <b>{</b>
<a name="223"><font color="#888888 "><i>  223   </i></font></a>          scale += h * f[i0 + 14 * i1] * c_b[i1];
<a name="224"><font color="#888888 "><i>  224   </i></font></a>        <b>}</b>
<a name="225"><font color="#888888 "><i>  225   </i></font></a>
<a name="226"><font color="#888888 "><i>  226   </i></font></a>        b_x[i0] = (x[i0] + scale) - sol2[i0];
<a name="227"><font color="#888888 "><i>  227   </i></font></a>      <b>}</b>
<a name="228"><font color="#888888 "><i>  228   </i></font></a>
<a name="229"><font color="#888888 "><i>  229   </i></font></a>      error_1 = 0.0;
<a name="230"><font color="#888888 "><i>  230   </i></font></a>      scale = 2.2250738585072014E-308;
<a name="231"><font color="#888888 "><i>  231   </i></font></a>      <font color="#112266"><b>for</b></font> (i = 0; i &lt; 14; i++) <b>{</b>
<a name="232"><font color="#888888 "><i>  232   </i></font></a>        absx = muDoubleScalarAbs(b_x[i]);
<a name="233"><font color="#888888 "><i>  233   </i></font></a>        <font color="#112266"><b>if</b></font> (absx &gt; scale) <b>{</b>
<a name="234"><font color="#888888 "><i>  234   </i></font></a>          b_t = scale / absx;
<a name="235"><font color="#888888 "><i>  235   </i></font></a>          error_1 = 1.0 + error_1 * b_t * b_t;
<a name="236"><font color="#888888 "><i>  236   </i></font></a>          scale = absx;
<a name="237"><font color="#888888 "><i>  237   </i></font></a>        <b>}</b> <font color="#112266"><b>else</b></font> <b>{</b>
<a name="238"><font color="#888888 "><i>  238   </i></font></a>          b_t = absx / scale;
<a name="239"><font color="#888888 "><i>  239   </i></font></a>          error_1 += b_t * b_t;
<a name="240"><font color="#888888 "><i>  240   </i></font></a>        <b>}</b>
<a name="241"><font color="#888888 "><i>  241   </i></font></a>      <b>}</b>
<a name="242"><font color="#888888 "><i>  242   </i></font></a>
<a name="243"><font color="#888888 "><i>  243   </i></font></a>      error_1 = scale * muDoubleScalarSqrt(error_1);
<a name="244"><font color="#888888 "><i>  244   </i></font></a>
<a name="245"><font color="#888888 "><i>  245   </i></font></a>      <font color="#117755"><i>/*  Estimate the error and the acceptable error */</i></font>
<a name="246"><font color="#888888 "><i>  246   </i></font></a>      scale = 0.0;
<a name="247"><font color="#888888 "><i>  247   </i></font></a>      i = 0;
<a name="248"><font color="#888888 "><i>  248   </i></font></a>      exitg2 = FALSE;
<a name="249"><font color="#888888 "><i>  249   </i></font></a>      <font color="#112266"><b>while</b></font> ((exitg2 == FALSE) &amp;&amp; (i &lt; 14)) <b>{</b>
<a name="250"><font color="#888888 "><i>  250   </i></font></a>        absx = muDoubleScalarAbs(x[i]);
<a name="251"><font color="#888888 "><i>  251   </i></font></a>        <font color="#112266"><b>if</b></font> (muDoubleScalarIsNaN(absx)) <b>{</b>
<a name="252"><font color="#888888 "><i>  252   </i></font></a>          scale = rtNaN;
<a name="253"><font color="#888888 "><i>  253   </i></font></a>          exitg2 = TRUE;
<a name="254"><font color="#888888 "><i>  254   </i></font></a>        <b>}</b> <font color="#112266"><b>else</b></font> <b>{</b>
<a name="255"><font color="#888888 "><i>  255   </i></font></a>          <font color="#112266"><b>if</b></font> (absx &gt; scale) <b>{</b>
<a name="256"><font color="#888888 "><i>  256   </i></font></a>            scale = absx;
<a name="257"><font color="#888888 "><i>  257   </i></font></a>          <b>}</b>
<a name="258"><font color="#888888 "><i>  258   </i></font></a>
<a name="259"><font color="#888888 "><i>  259   </i></font></a>          i++;
<a name="260"><font color="#888888 "><i>  260   </i></font></a>        <b>}</b>
<a name="261"><font color="#888888 "><i>  261   </i></font></a>      <b>}</b>
<a name="262"><font color="#888888 "><i>  262   </i></font></a>
<a name="263"><font color="#888888 "><i>  263   </i></font></a>      scale = auxdata-&gt;tol * muDoubleScalarMax(scale, 1.0);
<a name="264"><font color="#888888 "><i>  264   </i></font></a>
<a name="265"><font color="#888888 "><i>  265   </i></font></a>      <font color="#117755"><i>/*  Update the solution only if the error is acceptable */</i></font>
<a name="266"><font color="#888888 "><i>  266   </i></font></a>      <font color="#112266"><b>if</b></font> (error_1 &lt;= scale) <b>{</b>
<a name="267"><font color="#888888 "><i>  267   </i></font></a>        t += h;
<a name="268"><font color="#888888 "><i>  268   </i></font></a>        *tout = t;
<a name="269"><font color="#888888 "><i>  269   </i></font></a>        <font color="#112266"><b>for</b></font> (i = 0; i &lt; 14; i++) <b>{</b>
<a name="270"><font color="#888888 "><i>  270   </i></font></a>          x[i] = sol2[i];
<a name="271"><font color="#888888 "><i>  271   </i></font></a>          xout[i] = sol2[i];
<a name="272"><font color="#888888 "><i>  272   </i></font></a>        <b>}</b>
<a name="273"><font color="#888888 "><i>  273   </i></font></a>      <b>}</b>
<a name="274"><font color="#888888 "><i>  274   </i></font></a>
<a name="275"><font color="#888888 "><i>  275   </i></font></a>      <font color="#117755"><i>/*  Step control */</i></font>
<a name="276"><font color="#888888 "><i>  276   </i></font></a>      <font color="#112266"><b>if</b></font> (error_1 == 0.0) <b>{</b>
<a name="277"><font color="#888888 "><i>  277   </i></font></a>        error_1 = 2.2204460492503131E-15;
<a name="278"><font color="#888888 "><i>  278   </i></font></a>      <b>}</b>
<a name="279"><font color="#888888 "><i>  279   </i></font></a>
<a name="280"><font color="#888888 "><i>  280   </i></font></a>      emlrtPushRtStackR2012b(&amp;<a href="#var_c_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_c_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">c_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="281"><font color="#888888 "><i>  281   </i></font></a>      scale /= error_1;
<a name="282"><font color="#888888 "><i>  282   </i></font></a>      emlrtPushRtStackR2012b(&amp;<a href="#var_i_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_i_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">i_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="283"><font color="#888888 "><i>  283   </i></font></a>      emlrtPushRtStackR2012b(&amp;<a href="#var_j_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_j_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">j_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="284"><font color="#888888 "><i>  284   </i></font></a>      <font color="#112266"><b>if</b></font> (scale &lt; 0.0) <b>{</b>
<a name="285"><font color="#888888 "><i>  285   </i></font></a>        emlrtPushRtStackR2012b(&amp;<a href="#var_k_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_k_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">k_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="286"><font color="#888888 "><i>  286   </i></font></a>        <a href="#fcn_b_eml_error" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_b_eml_error');" target="rtwreport_document_frame"><font color="#1122aa">b_eml_error</font></a>();
<a name="287"><font color="#888888 "><i>  287   </i></font></a>        emlrtPopRtStackR2012b(&amp;<a href="#var_k_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_k_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">k_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="288"><font color="#888888 "><i>  288   </i></font></a>      <b>}</b>
<a name="289"><font color="#888888 "><i>  289   </i></font></a>
<a name="290"><font color="#888888 "><i>  290   </i></font></a>      emlrtPopRtStackR2012b(&amp;<a href="#var_j_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_j_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">j_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="291"><font color="#888888 "><i>  291   </i></font></a>      emlrtPopRtStackR2012b(&amp;<a href="#var_i_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_i_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">i_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="292"><font color="#888888 "><i>  292   </i></font></a>      h = muDoubleScalarMin(hmax, 0.9 * h * muDoubleScalarPower(scale, 0.125));
<a name="293"><font color="#888888 "><i>  293   </i></font></a>      emlrtPopRtStackR2012b(&amp;<a href="#var_c_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_c_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">c_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="294"><font color="#888888 "><i>  294   </i></font></a>      emlrtBreakCheckFastR2012b(<a href="ode78__data_c.html#var_emlrtBreakCheckR2012bFlagVar" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtBreakCheckR2012bFlagVar');" target="rtwreport_document_frame"><font color="#1122aa">emlrtBreakCheckR2012bFlagVar</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="295"><font color="#888888 "><i>  295   </i></font></a>      <font color="#112266"><b>if</b></font> (muDoubleScalarAbs(h) &lt;= 2.2204460492503131E-16) <b>{</b>
<a name="296"><font color="#888888 "><i>  296   </i></font></a>        exitg1 = TRUE;
<a name="297"><font color="#888888 "><i>  297   </i></font></a>      <b>}</b>
<a name="298"><font color="#888888 "><i>  298   </i></font></a>    <b>}</b>
<a name="299"><font color="#888888 "><i>  299   </i></font></a>  <b>}</b>
<a name="300"><font color="#888888 "><i>  300   </i></font></a><b>}</b>
<a name="301"><font color="#888888 "><i>  301   </i></font></a>
<a name="302"><font color="#888888 "><i>  302   </i></font></a><font color="#117755"><i>/* End of code generation (ode78_.c) */</i></font>
<a name="303"><font color="#888888 "><i>  303   </i></font></a></pre>
</td></tr></table>
</p>
</body>
</html>